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© Signal to density mapping with controlled contrast. 

© A method of determining a conversion of a digi- 
tal image representation of a radiologic image into 
density values wherein the gradation curve of a 
conventional screen/film system is approximated. 
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Field of the invention. 

The present invention is in the field of digital 
image processing. The invention more specifically 
relates to a method of creating a gradation curve, 
i.e. a curve representing the relationship between 
the optical density of a reproduced image and the 
corresponding digital signal value. 

Background of the invention 

In the field of digital radiography a system has 
been developed wherein x-rays transmitted by an 
exposed object are stored on a photostimulable 
phosphor screen. 

The photostimulable phosphor screen is read- 
out by stimulating it with laser light so that it 
liberates its latent image by the emission of image- 
wise modulated light. The emitted light is then 
detected and converted into an analog electrical 
representation for example by means of a 
photomultiplier. The signal is then digitized, pro- 
cessed and subsequently used to control the hard 
copy recording in image reproducing system such 
as a laser recorder. For diagnostic purposes the 
reproduction is viewed and analysed on an light- 
box. 

Examples of photostimulable phosphors that 
can be used are extensively described in European 
patent application 91200511.3 filed on 8.03.91. 

In conventional radiographic systems the op- 
tical density on the film is a function of the expo- 
sure dose. The variation of the density relative to 
the exposure dose (or the logarithm of the expo- 
sure dose), commonly referred to as the gradation 
curve or the sensitometric curve of the recording 
material is for each type of recording material 
optimized with respect to the application for which 
that type of material was developed. 

Optimization of the gradation curve i.a. means 
that the slope of said curve, that is a measure of 
the contrast rendition in each density region, is 
high enough so as to enable discerning the most 
subtle detail when examining the reproduction on a 
light-box taking into account the perceptibility char- 
acteristics of the human eye as a function of the 
intensity, the available density range and the rela- 
tive diagnostic importance of different densities in 
the image to be reproduced. 

This aspect is different in digital radiography 
since in such a system the recording medium and 
reproduction medium are not identical. Therefore a 
relation has to be explicitly defined between the 
digital signal representative of the light emitted by 
the phosphor upon radiation and the corresponding 
density on the recording material. This relation will 
furtheron be referred to as the digital gradation 
curve. 



Implementation of a digital gradation curve is a 
relatively simple task. Due to the large number of 
degrees of freedom and the large number of pa- 
rameters that can be set many gradation curves 
5 suited for different types of applications can be 
defined. 

However, when such a large number of grada- 
tion curves needs to be stored in an electronic 
memory a lot of storage capacity is required. 

io In European patent specification EP 0 110 185 

this problem has been described. According to this 
disclosure the problem is solved by a method of 
obtaining many desired gradation curves from 
about one of several (about ten) basic gradation 

75 curves by rotating or parallel shifting the basic 
gradation correction curve according to the object, 
the radiation source voltage, the purpose or the 
like. 

US patent 4,276,473 issued June 30, 1981 also 

20 relates to the conversion of digital signal values 
representing a radiologic image into corresponding 
density values. The level of the electric signal is 
converted in such a way that the maximum value of 
the electric signal corresponding to maximum den- 

25 sity of the radiation image is converted to a level 
resulting in a certain value of the optical density of 
the reproduced image on the recording material. 
The minimum level corresponding to the minimum 
density is converted to a level resulting in the fog 

30 density of the recording material. 

Hence a signal range that is entirely dependent 
on the read-out image information is mapped onto 
a fixed density range. When however the latitude of 
the read -out signal range varies, for example be- 

35 cause the fysionomy (the thickness) of the exposed 
patients varies, the average contrast in the re- 
productions will also vary. There is no control on 
the variation of the contrast obtained in hard copies 
of images with different information content. 

40 

Objects of the invention 

It is an object of the present invention to pro- 
vide a method for image reproduction wherein a 

45 gradation curve is used for converting digital image 
signals obtained by reading out an image stored in 
a stimulable phosphor screen into corresponding 
density values in the reproduction. 

It is a further object to provide such a method 

so wherein said curve is defined by means of a set of 
parameters that are familiar to the radiologist who 
is acquainted with classical radiology wherein the 
recording material and the reproduction material 
are the same conventional photographic film for x- 

55 ray applications. 
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Statement of the invention 

In accordance with the present invention the 
above objects are achieved by an image reproduc- 
ing method comprising the steps of scanning a 
stimulable phosphor screen storing a radiation im- 
age with stimulating rays, detecting the light emit- 
ted upon stimulation and converting it into a digital 
signal representation and creating a gradation 
curve for mapping the signal values of said digital 
signal representation onto corresponding density 
values and recording a hard copy image on a 
recording material taking into account said grada- 
tion curve, characterised in that said gradation 
curve is created by the steps of 

- defining a canonical monotonical non-con- 
stant function Y(x) in an orthogonal coordi- 
nate system (x,y) between arbitrary values xo 
and xt and y min and y max , 

- determining minimum and maximum density 
values D sm m and D smox , 

- determining minimum and maximum signal 
values S min and S max constituting a range 
wherein the gradation is defined by the ca- 
nonical form and outside of which signal val- 
ues will be mapped onto D smin or D smax , 

- defining the gradation curve as 

D(s) = D smin + (Y(x(s)) - y min ) . (D smax - 
) / (Ymax - y m in) wherein D(s) is the den- 
sity in the reproduction as a function of the 
signal value s, 

x(s) = xo if s is smaller than or equal to S min , 
x(s) = xo + (s-S mjn ).(xi -xo )/(S max -S mjn ) if s is 
comprised between S min and S max , and 
x(s) = xi if s is larger than or equal to S max - 
A gradation curve is a curve representing the 
conversion of digital signal values into correspond- 
ing density values envisioned in the hard copy. 

In a preferred embodiment the parameters 
Smin, Smax, D smax and tne canonical function Y(x) 
are determined as a function of the type of exami- 
nation performed. D smin is determined by the char- 
acteristics of the recording material used. In the 
case of a photographic recording material it is 
about equal to the fog density (e.g. D = 0.2). 

The shape of the canonical function Y(x) pref- 
erably approximates the shape of a gradation curve 
of a conventional photographic x-ray film the ra- 
diologist is used to work with. 

When the definition of a specific gradation 
curve representing a conversion from digital signal 
values into density values is performed starting 
from a predetermined function in a canonical form, 
then a calibration of the two axis of an orthogonal 
coordinate system wherein the canonical form is 
defined, is required, (cfr. Fig. 3 referred to 
herinafter) 



The calibration of the density axis is deter- 
mined by the values D smln and D smax that are a 
function of the inherent characteristics of the re- 
cording material and of the type of examination. 

5 The signal axis of the digital gradation curve 

that is defined is calibrated by S min and S ma x- In 
between these values the conversion of digital sig- 
nal values into density values is determined by the 
canonical function, signal values smaller than S mi n 

10 will be mapped onto D smin and signal values larger 
than S max will be mapped onto D smax . This trans- 
formation is obtained by means of the function x(s) 
as defined hereinbefore. 

It was the inventor's object to provide a digital 

ts gradation curve that is characterised by parameters 
that are familiar to the radiologist who is exper- 
ienced in the field of conventional radiography. 

A specific photographic quantity that is very 
familiar to the radiologist is the "exposure latitude" 

20 of a conventional radiographic film. A typical lati- 
tude value is 1 .5. 

Therefore in a preferred embodiment of a 
method as hereinbefore described the range be- 
tween S mjn and S max is constant and equal to a 

25 predetermined envisioned latitude L in the hard 
copy, so S max = S mi n + L 

Once the range between S m m and Smax is de- 
termined, it is still required to determine the posi- 
tion of this determined range on the signal axis, 

30 more specifically the position of this range relative 
to the range of the image signal values obtained by 
reading out the photostimulable phosphor screen. 

According to this invention this is performed by 
taking into account the statistical information of the 

35 read-out image signal, more specifically the image 
histogram. 

By evaluation of the image histogram the rel- 
evant signal range being the range comprising in- 
formation that is relevant for diagnostic purposes is 

40 determined. 

Then an alignment is performed of a signal 
value located at a predetermined relative position 
within the range defined by S mjn and S max wherein 
the canonical function is applicable with a point 

45 located at the same predetermined relative position 
within the relevant signal range. For example 62 % 
of the signal range whereto the conversion via the 
canonic function is confined is aligned with 62 % of 
the relevant signal range. 

so For reproducing images represented by a digi- 

tal signal representation the mapping of a fixed 
predetermined range onto corresponding density 
values is advantageous relative to the mapping of a 
data-dependent range because in the former case 

55 the gradation obtained in the hard copy will be 
constant. The information that is important for the 
radiologist when making a diagnosis is contained in 
the image contrast (i.e. in density differences), so 
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the control of the contrast rendition is most impor- 
tant This direct control of the contrast rendition is 
obtained by the method of the present invention. 

In the alternative case in which the range that 
is mapped onto density values is without any limi- 
tations dependent on the image information, a 
change of the signal range brings about an implicit 
rescaling of the gradation curve. This has an effect 
on the average contrast in the hard copy. As a 
result hereof the gradation of the hard copy image 
may for example vary considerably in dependence 
on the thickness of the irradiated patient. 

However, the use of a fixed signal range has 
disadvantages too. Indeed, when the range of di- 
agnostically relevant data determined by evaluation 
of the image histogram is larger than the fixed 
signal range to which the conversion is confined, a 
part of the diagnostically relevant data will be 
mapped onto the minimum or maximum density so 
that the local image contrast in these parts of the 
image is equal to zero. 

To solve this problem the present invention 
provides an alternative embodiment of the de- 
scribed method wherein the range between Sm in 
and S max , earlier referred to as the latitude, that is 
mapped onto density values in between D sm i n and 
D smax is adjusted to the dimensions of the range 
wherein relevant histogram data are found within 
given margins so that the corresponding variations 
of the slope of the gradation curve determined as 
hereinbefore described, will be controlled. 

The choice between the two above described 
particular embodiments of the present invention is 
preferably made in dependence of the type of 
examination and will commonly be provided for the 
user in the form of a menu on screen. 

A hard copy image can be obtained by scan- 
ning a photographic recording material with a laser 
beam that is modulated so as to produce in each 
pixel on the recording material the specific density 
value that corresponds with the read-out signal 
value in said pixel taking into account the above- 
described gradation curve. 

Other recording methods such as thermal sub- 
limation transfer recording may be envisioned. The 
process of thermal sublimation transfer comprises 
placing a dye layer of a donor element in face-to- 
face relation with a dye receiving layer of a re- 
ceiver sheet and image-wise heating from the back 
of the donor element. 

For supplying heat energy thermal heads, laser 
light, infrared flash or heated pens can be used. 

Alternatively the support of the dye-donor ele- 
ment may be an electrically resistive ribbon 
wherein current is injected by electrically address- 
ing a print head electrode resulting in a highly 
localised heating of the ribbon beneath the relevant 
electrode. 



Brief description of the drawings 

Particular aspects of the present invention as 
well as preferred embodiments thereof will be ex- 
5 plained by means of the corresponding drawings 
wherein 

Fig. 1 generally shows a system for reading out 
an image stored in a photostimulable phosphor 
screen, 

10 Fig. 2 is a block sheme illustrating the data flow, 
Fig. 3 illustrates the method of defining a grada- 
tion curve starting from a given canonical func- 
tion Y(x). 

Fig. 4a-b illustrate the alignment of the range 
75 confined by S mln and S max with the diagnosti- 
cally relevant signal range defined by So and 
S,. 

Fig. 5 and Fig. 6 are examples of canonical 
functions that can be used for creating a grada- 

20 tion curve. 

Figure 1 generally shows an apparatus in 
which the method of the invention can be applied. 

A radiation image of an object was recorded on 
a photostimulable phosphor screen (3) by exposing 

25 (2) said screen to x-rays transmitted through the 
object (not shown). The stimulable phosphor 
screen was conveyed in a cassette (3) provided 
with an electrically erasable programmable read 
only memory (EEPROM). In an identification station 

30 4 various kinds of data, for example patient iden- 
tification data (name, date of birth) and data relat- 
ing to the exposure and/or to the signal processing 
were recorded onto the EEPROM. 

In a radiation image read-out apparatus 1 the 

35 latent image stored in the photostimulable phos- 
phor screen was read-out by scanning the phos- 
phor screen with stimulating rays emitted by a 
laser. The stimulating rays were deflected into the 
main scanning direction by means of galvanometric 

40 deflection. The subscanning was performed by 
transporting the phosphor screen in the subscan- 
ning direction. The stimulated emission was di- 
rected onto a photomultiplier for conversion into an 
electrical image representation. 

45 The subsequent data flow is illustrated in figure 

2. The output signal of the photomultiplier was 
converted into a logaritmic quantity log E (E being 
the exposure value), and next the signal was quan- 
tised. This quantised image signal, called the raw 

so image signal, was sent to the image processor 
(figure 1, numeral 7) where it was stored in an 
internal buffer. Without any modifications it was 
also sent from the image processor to the review 
console 5 where it was temporarily stored on a 

55 hard disc. This back up ensured that the signal was 
never lost, even not when any of the components 
of the apparatus would fail and that the signal could 
be retrieved for any kind of later processing, for 
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example processing with different parameter set- 
ting. This feature could be used when the result of 
the on-line processing was unsatisfactory due to 
bad exposure conditions or inadequate selection of 
the processing parameters. 

Next the curve representing the mapping of 
signal values (in this case log exposure) values into 
corresponding density values was to be deter- 
mined. 

The way in which this gradation curve was 
defined in accordance with the present invention is 
illustrated by means of figure 3 wherein the quad- 
rants are numerated in clockwise order starting 
from the upper left corner. 

First the minimum density value D smln and the 
maximum density value D smax envisioned in the 
hard copy were defined, D smjn was equal to fog 
density and D smax was equal to 3.0. 

Then a canonical function defined in an or- 
thogonal coordinate system between x Q . xi and 
Vmin.ymax was retrieved from the internal memory of 
the signal processor. An example of such a curve 
is depicted in the fourth quadrant of figure 3. Spe- 
cific embodiments are shown in figures 5 and 6. 

Next two values S min and S max were deter- 
mined that constitute a range wherein the conver- 
sion of signal values onto density values is deter- 
mined by the specific shape of the canonical func- 
tion. Signal values smaller than S min are mapped 
onto D smjn , signal values greater than D smax are 
mapped onto D smax . 

In this embodiment the latitude of said range 
was a fixed value L = 1.5 log exposure 
(corresponding with the latitude of a conventional x- 
ray film the radiologist is used to work with) and 
the position of S mln was determined relative to the 
diagnostically relevant signal range. S max was then 
calculated as S mjn + L. 

For determining the position of S min relative to 
the relevant signal range, the extreme values So 
and Si of the diagnostically relevant signal range 
were first determined by evaluation of the image 
histogram. 

Then a small offset dSi = 0.3 log E was added 
to Si. This ensures that the density in the hard 
copy corresponding with the maximum value of the 
diagnostically relevant signal range does not be- 
come too dark. 

The positioning of the range S max -S mm relative 
to the range Si -So was performed by aligning a 
fraction of the latter range with the same fraction of 
the former range. 

Mathematically this fraction can be expressed 
as A(Si +dSi-So-dSo). 

Then the alignment can be formulated math- 
ematically as S min = So + dSo+ A(Si +dSi-So-dSo) 
- A.L and S max = S m i n + L J A being an integer value 
greater than or equal to 0 and smaller than or equal 



to 1. 

This mathematical formulation is illustrated in 
figure 4a for A = 0.62. 

The embodiment wherein A was equal to 1 and 
5 dSo was equal to zero is illustrated in figure 4b. 

Next the curve is calculated as D(s) = D smin + 
(Y(x(s)) - y mln ). (D smax - D smjn ) / (y max - y mjn ) wherein 
D(s) is the density in the reproduction as a function 
of the signal value s and wherein x(s) = xo if s is 
to smaller or equal than S mjn , x(s) = xo + (s-S m in)(xi - 
xo)/(S max -S min ) if s is comprised between S mjn and 
S m ax. and x(s)=xi if s is larger or equal to S max . 
The function x(s) is shown in figure 3, second 
quadrant. 

is Finally the signal is processed on the basis of 

this gradation curve and the processed signal is 
then applied to a digital to analog converter, the 
output of the converter modulates the output of a 
laser recorder 6 for making a hard-copy on photo- 

20 graphic film. 

This kind of gradation processing has been 
applied to a series of radiation images of chest 
radiographs. 

The alignment of the range S min to S max relative 

25 to the diagnostically relevant range between So and 
Si is performed as hereinbefore described. In 
these reproductions a quasi constant density was 
obtained in the darkest part of the lungs and the 
change of contrast in the remaining part of the 

30 image is identical in each of the reproductions, it 
was independent of the thickness of the patient that 
was irradiated. The canonical function shown in 
figure 5 was used. The shape of this function 
approximates the gradation curve of a known pho- 

35 tographic film and the outlook of the reproduction 
was very familiar to the radiologist. 

The addition of a small offset value dSi to Si 
as described hereinbefore provides that the density 
in the darkest part of the lungs does not become 

40 too high. 

It was also possible to subtract a small offset 
value dSo from So, in this example dSo was equal 
to zero. 

Alternatively the latitude between S m j n and S max 
45 was not fixed but was adapted to the latitude of the 
diagnostically relevant signal range determined by 
evaluation of the histogram. The adaption was how- 
ever limited within predetermined margins defined 
as a minimum latitude L min and a maximum latitude 
50 L max . 

Mathematically this can be formulated as fol- 
lows: 

Smax = S mjn + L and S min = So + dS 0 + A.(Si+dSi-S 
55 o-dSo) - A.L 

wherein 

L = L min if L mi n>Si +dSi-So-dS 0 ; 
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L = Si + dSi -So -dSo if L min < = Si + dSi -So - 
dSo< = L max ; 

L = L max if Si +dSi-So-dSo>S max - 

and wherein -L max /2< =dSo< = L max /2 and 

L ma x/2< = dSlo< = L max /2. 

This embodiment was applied for the reproduc- 
tion of a series of chest radiographs. The set 
values for the different parameters were the follow- 
ing: 

dS o =0; dSi=0.1 log E, A = 0, L min = 1 .2 log E; 
Lmax = 1-8 log E; D smax = 3.0; the canonical function 
that was used is shown in figure 6. 

In the reproductions of chest radiographs a 
quasi constant density was obtained in the me- 
diastinum (D = 0.2). The contrast in the mediasti- 
num and in the lungs is good in every reproduction 
provided that the effective latitude between the 
minimum and the maximum of the diagnostically 
relevant data determined by evaluation of the histo- 
gram is smaller than 1.8. In case this effective 
latitude exceeds 1.8, the contrast in the darkest 
parts of the lungs becomes zero. 

Hence, whether the embodiment of the present 
invention characterised by a fixed latitude or the 
embodiment characterised by a limited-variable 
latitude is selected depends on the subject wherein 
the radiologist is interested. In the present applica- 
tion one should preferably choose the embodiment 
incorporating a fixed latitude when the field of 
interest is situated in the lungs, whereas the em- 
bodiment incorporating a limited-variable latitude is 
to be chosen when the field of interest is situated 
in the mediastinum. 

Claims 

1. An image reproducing method comprising the 
steps of scanning a stimulable phosphor 
screen storing a radiation image with stimulat- 
ing rays, detecting the light emitted upon stim- 
ulation and converting it into a digital signal 
representation and creating a gradation curve 
for mapping the signal values of said digital 
signal representation onto corresponding den- 
sity values and recording a hard copy image 
on a recording material taking into account 
said gradation curve, characterised in that said 
gradation curve is created by the steps of 

- defining a canonical monotonical non- 
constant function Y(x) in an orthogonal 
coordinate system (x,y) between arbitrary 
values xo and xi and y min and y max , 

- determining minimum and maximum 
density values D smin and D smax , 

- determining minimum and maximum sig- 
nal values S m in and S max constituting a 
range wherein the gradation is defined 
by the canonical form and outside of 



which signal values will be mapped onto 

Dsmin Or D smax , 

- defining the gradation curve as 

D(S) = D smin + (Y(X(S)) - y m in) - (Dsmax " 

5 D smln ) / (y max - y min ) wherein D(s) is the 

density in the reproduction as a function 
of the signal value s, 

x(s) = x 0 if s is smaller than or equal to 

Snnin> 

10 X(S) = X 0 +(S-S m i n ).(X 1 -Xo)/(S m ax-S mi n) if S 

is comprised between S min and S max » and 
x(s) = xi if s is larger than or equal to 

Smax- 

75 2. A method according to claim 1 wherein Y(x), 
S mjn , S max and D smax are determined as a func- 
tion of an examination type and D smjn is equal 
to the fog level of the recording material. 

20 3. A method according to claim 1 wherein 
S max = S mln + L, L is a constant latitude value 
envisioned in a reproduction of said image and 
S mjn is determined by aligning a signal value at 
a predetermined relative position between S min 

25 and S maix with a signal value at the same 

predetermined relative position between So- 
dS 0 and Si +dSi wherein S1-S0 is the di- 
agnostically relevant signal range and wherein 
dSo and dSi are offset values comprised in 

30 the range confined by -L/2 and L/2. 

4. A method according to claim 1 wherein the 
signal is converted into a logarithmic quantity 
log exposure before being digitized and 

35 wherein dSo and dSi are situated between -0=3 

log exposure and 0.3 log exposure. 

5. A method according to claim 1 wherein 
S max = S min + L and S mi n = So + dS 0 + A(Si + dSi - 

40 S 0 -dSo)-A.L 

wherein L is equal to a predetermined value 
Lmin if U, in is greater than Si + dSi-S 0 -dSo or 
equal to a predetermined value L ma x if Lmax is 
smaller than Si +dSi-S 0 -dS 0 or L is equal to 

45 Si +dSi-S 0 -dSo if Si +dSi-So-dS 0 is com- 

prised between L min and L maK and wherein A is 
a constant value between 0 and 1, dSo and 
dSi are offset values comprised between - 
Lmax/2 and L max /2 

50 

6. A method according to claim 5 wherein Y(x), 
Smin, Smax and D smax are determined as a func- 
tion of an examination type, and D smm is equal 
to the fog level of the recording material. 

55 

7. A method according to claim 5 wherein the 
signal is converted into a logarithmic quantity 
log exposure before being digitized and 
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wherein dS 0 and dSi are situated between -0.3 
log exposure and 0.3 log exposure. 
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